arXiv:1502.07071vl [quant-ph] 25 Feb 2015 


Observation of a phononic Mollow triplet in a hybrid spin-nanomechanical system 
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Reminiscent of the bound character of a qubit’s dynamics confined on the Bloch 
sphere, the observation of a Mollow triplet in the resonantly driven qubit fiuorescence 
spectrum represents one of the founding signatures of Quantum Electrodynamics. Here 
we report on its observation in a hybrid spin-nanomechanical system, where a Nitro¬ 
gen Vacancy spin qubit is magnetically coupled to the vibrations of a Silicon Carbide 
nanowire. A resonant microwave field turns the originally parametric hybrid interac¬ 
tion into a resonant process, where acoustic phonons are now able to induce transitions 
between the dressed qubit states, leading to synchronized spin-oscillator dynamics. We 
further explore the vectorial character of the hybrid coupling to the bidimensional de¬ 
formations of the nanowire. The demonstrated microwave assisted synchronization 
of the spin-oscillator dynamics opens novel perspectives for the exploration of spin- 
dependent forces, the key-ingredient for quantum state transfer. 


Introduction- A mechanical oscillator coupled to a two 
level system (TLS) is a versatile basis to study the in¬ 
teraction between macroscopic and purely quantum ob¬ 
jects. This unconventional combination mm is a promis¬ 
ing route towards the generation of non-classical states 
of motion of macroscopic objects. Hybrid coupling signa¬ 
tures have now been demonstrated between a mechanical 
oscillator and Bose Einstein Condensates 0 g , super¬ 
conducting qubits m, solid state single spins EHni, 
molecules [18] or quantum dots [T9H23] . 

The hybrid interaction coupling phonons and qubits is in 
profound analogy with quantum electrodynamics (QED) 
where hallmark experiments revealing the interplay be¬ 
tween atoms and photons have permitted exploring the 
foundations of quantum mechanics. In particular, the ap¬ 
parition of a Mollow triplet in atomic fiuorescence spectra 
[25] . characterized by the onset of sidebands appearing 
on each side of the pump frequency with splitting propor¬ 
tional to the laser field strength, is one of the character¬ 
istic signatures of the strongly driven Jaynes-Gumming 
interaction. Along with the Autler-Townes doublet [26| 
or vacuum Rabi oscillations, it expresses the dressing of 
the atom with the optical photon field m- Since then, 
Mollow triplets were observed in atomic vapors [28l |29] , 
single molecules [301EB, single quantum dots [32j |33] or 
superconducting qubits [34] coupled to photon fields in 
the optical or microwave (MW) domain. Here we report 
on the observation of a phononic Mollow triplet, where 
the phonon field of a nanomechanical oscillator dresses a 
MW dressed single spin qubit. 

Our hybrid device consists of a single Nitrogen Vacancy 
(NV) spin qubit hosted in a diamond nanocrystal at¬ 
tached to the vibrating extremity of a Silicon Garbide 
(SiG) nanowire [TO]. A strong magnetic field gradi¬ 
ent couples both components through a spatially depen¬ 
dent Zeeman effect (see Eig. la). Eormally, the generic 
hybrid spin-oscillator Hamiltonian can be expressed as 
^)Fhgi{aFa^)o-i, where ujqI2ii is the 


qubit energy, FL^/2ii the oscillator frequency, a^{a) the 
phonon creation (annihilation) operator, (Ji the Pauli ma¬ 
trices of the spin qubit quantized along the 2 ; axis and Qi 
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FIG. 1. A hybrid spin qubit-nanomechanical system. 

a. A single NV spin qubit hosted in a diamond nanocrystal 
is attached at the oscillating extremity of a suspended SiC 
nanowire. A strong magnetic field gradient source is micro- 
positioned in the vicinity of the hybrid system to magnetically 
couple the spin state to the vibrations of the nanoresonator 
through the Zeeman effect, b. In the parametric coupling 
regime, the mechanical motion modulates the qubit energy 
with an amplitude (5ct;o. c. Scanning fluorescence image of the 
hybrid device, d. Autocorrelation function of the NV spin 
qubit fluorescence revealing the presence of a single defect, 
e. ESR spectroscopy of the suspended spin qubit in a weak 
magnetic field, highlighting the spin state dependence of the 
average emitted fluorescence (FWHM~ 4MHz), e. Brownian 
motion of the nanowire measured in reflection (below) reveal¬ 
ing the two fundamental eigenmodes, which can be coherently 
driven through electrostatic or piezo actuations. The response 
curve (above) permits to determine the local orientation gf 
of the force vector as well as its magnitude EH. 
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the respective coupling constants. As in quantum electro¬ 
dynamics, a distinction between the resonant and para¬ 
metric regimes can be made. In the first situation, de¬ 
scribed by the interaction Hamiltonian hgx^y{ci-\-a^)(Jx^y^ 
the mechanical oscillator and qubit can coherently ex¬ 
change single excitations if they have similar frequen¬ 
cies. Several hybrid mechanical systems are exploring 
this regime, either through a direct interaction P l35l[36] 
or mediated by a bus cavity 13 EH Eg. In the case of 
parametric hybrid coupling, the mechanical motion mod¬ 
ulates the qubit energy according to the coupling Hamil¬ 
tonian hgz{a + o))(Jz' Reciprocally, the qubit exerts a 
state dependent force on the oscillator which represents 
the key ingredient for quantum state transfer between 
both components. This configuration enables hybrid cou¬ 
pling between components with totally different excita¬ 
tion energies. Moreover by employing a resonant MW 
tone driving Rabi precession of the qubit at frequencies 
flR/27r close to the mechanical frequency, it is possible 
to let both components of the hybrid system evolve on 
similar time scales [15] [39] . By doing so the parametric 
interaction with the original qubit is turned into a reso¬ 
nant coupling to the MW dressed qubit. This configura¬ 
tion enables the observation of a phononic Mollow triplet, 
provided that the phonon field is coherently driven and 
that the oscillator frequency is larger than the qubit de¬ 
cay rate Tgpin, which corresponds to the so called resolved 
sideband regime > Tspin) [IS] of the parametric in¬ 
teraction. In that situation we introduce the dynamical 
parametric modulation strength Sujq^ which denotes the 
classical amplitude of the mechanically driven parametric 
modulation (see Fig. lb). 

The experimental setup is sketched in Fig. la. The 
nanomechanical oscillator is a 6 jam long SiC nanowire of 
300 nm diameter, suspended from a sharp metallic tip. 
Its moving extremity is functionalized with a 50 nm 
nano-diamond hosting a single NV defect. The hybrid 
system is investigated with a confocal microscope appa¬ 
ratus (see SI) and a 532 nm pump laser. It serves for both 
measuring the vibrations of the nanowire using either the 
transmitted or reflected light beams and for optical po¬ 
larization and readout of the qubit using spin state de¬ 
pendent fluorescence detection [40] . 

Measurements of the nanowire Brownian motion permit 
determining the mechanical properties of the fundamen¬ 
tal flexural eigenmodes [24]. These are aligned along 
two perpendicular directions ei, 2 , see Fig. If, tilted by 
^ 10° with respect to the optical axis, at frequencies 
^m/27r (m=l or 2) of 5.99 and 6.29 MHz respectively, 
with mechanical damping rates ri(' 2 )/ 27 r = 180 (190) kHz 
in air, limited by acoustic emission. The measured ef¬ 
fective masses of Meff ~ 10“^^ kg correspond to a spa¬ 
tial spreading of their Brownian motion over Axth = 
~ 52 pm with zero-point fluctuations 
of Axq ~ 36 fm. Using a resonant force actuation (5F, 
either piezoelectric or electrostatic, it is possible to drive 


vibrations of the the nanowire around its rest position 
tq. Its vectorial deflection (5r(t) can thus be expressed in 
frequency space as: 

Sr[n] = XmM {SF,^ + (5F[n] • e^) (1) 

m=l,2 

using the mechanical susceptibilities l/xrn[^] = 
Meff (Um^ — — iUFm) and the independent Langevin 

forces SF^. By adjusting the drive frequency as well as 
the orientation ep of the force vector with respect to the 
nanowire eigenmodes it is thus possible to generate differ¬ 
ent trajectories in the oscillation plane. This will permit 
exploring the vectorial character of the magnetic hybrid 
interaction. 

The backscattered fluorescence of the NV defect is col¬ 
lected by avalanche photodiodes arranged in a Hanbury 
Brown and Twiss configuration which allows to confirm 
the presence of a single NV defect through fluorescence 
intensity autocorrelation functions (see Fig. Id.) m- 
The qubit can be initialized in its ground state through 
optical pumping, manipulated with quasi-resonant MW 
fields and readout exploiting its spin state dependent 
fluorescence rate uni- When immersing the suspended 
NV defect in a strong magnetic field gradient the spin 
state becomes position dependent, as described by the 
Zeeman Hamiltonian —ggs^ • B(ro + ^r). At first or¬ 
der the nanomechanical oscillator’s vectorial deflection 
and the spin are dynamically coupled through a magnetic 
field gradient according to the interaction Hamiltonian: 

i^int =-5MB cr ■ (<^r ■ V)B|r^ . (2) 

Thus, it is generally possible to tune the interaction from 
resonant {<Jx,y) to parametric (az) by adjusting the qubit 
frequency and the topography of the magnetic field gra¬ 
dient. In the following we will only investigate the para¬ 
metric case, where the nanomotion modulates the qubit 
energy. 

A 18 gm diameter NdFeB hard magnetic sphere gen¬ 
erates the magnetic field gradient needed to generate a 
large hybrid coupling strength. It is positioned onto a 
narrow gold stripline antenna, see Fig. 2b, used to deliver 
the resonant microwave field and especially designed to 
fit between the high numerical aperture objectives with¬ 
out compromising the photon collection efficiency (see SI) 
and the spatial access to the magnetic gradient source. 
The latter can be piezo scanned with respect to the sus¬ 
pended NV spin with nanometric precision. The follow¬ 
ing restrictions dictate the necessary configuration of the 
experimental setup. Rotating the nanowire while opti¬ 
mizing the collected fluorescence allows to align the in¬ 
trinsic spin’s quantization direction along the optical z 
axis m- The remanent magnetization of the NdFeB 
particle was oriented in a strong magnetic field (~ 1.5 T) 
and aligned with the NV natural quantization axis to 
preserve the spin-selective readout efficiency. The mag- 
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netic bead stray field also allows to polarize the nu¬ 
clear spin by working at the excited state level anticross¬ 
ing (B ^ 50niTez) [42l HS], which permits to restrict 
our system to a pure TLS through MW frequency selec¬ 
tion. Finally, the last fundamental requirement consists 
in reaching a large parametric coupling strength, orient¬ 
ing 'VBz with the direction of the eigenmode of interest. 

To determine the vectorial parametric coupling 
strength, the spatial dependence of the spin qubit en¬ 
ergy u;o(r) was measured by collecting the NV fluores¬ 
cence while scanning the position of the micromagnet in 
presence of a continuous microwave tone. Typical fluo¬ 
rescence maps are shown in Fig. 2d for varying MW fre¬ 
quencies. The projections of the qubit iso-energy sur¬ 
faces on the oscillation plane, see Fig. 2d, appear as dark 
resonant slices [44f[T7] . In addition, a global fluores¬ 
cence quenching which indicates regions with a strong 



FIG. 2. Mapping the parametric vectorial hybrid 
interaction, a. The coupling strength is determined by 
measuring the spatial dependence of the spin qubit energy 
huJo(r) in absence of mechanical drive for different positions 
of the micromagnet, b. Colored SEM image showing the 
30 /xm wide suspended MW waveguide supporting the NdFeB 
microsphere, which is piezo-scanned in the vicinity of the NV. 
c. Fluorescence rates measured at different horizontal posi¬ 
tions (22, 20 and 17.5 /xm) while scanning the micromagnet in 
the vertical xy plane. The mechanism of fluorescence quench¬ 
ing in transverse magnetic field strongly confines the suitable 
working points, d. Similar measurements in the horizontal 
xz plane in presence of a single MW tone at varying frequen¬ 
cies. e. Corresponding numerical simulations based on a 
purely dipolar magnetic field (see SI). The resonant slices ap¬ 
pear as dark lines which identify the iso-spin energy surfaces 
h(jUo(r), whose spatial derivative permits determining the vec¬ 
torial parametric coupling strength, A = Vcjq. Its projection 
onto the oscillation plane is shown in f. as a color map while 
its orientation is indicated by the orange arrows. The dashed 
lines reproduce the iso-energy lines (100 MHz spacing). O 
indicates the working point used in the following. 


off-axis magnetic held [48]. Reproducing this measure¬ 
ment for varying microwave frequencies permits deter¬ 
mination of co’o(r) (see Fig. 2f). When moving in those 
strong magnetic held gradients, the suspended spin qubit 
undergoes a dynamical parametric energy modulation of 
SuJo{t) = Sr{t) • VciJoIro? which is determined by evaluat¬ 
ing the gradient of the iso-energy map in the ( 01 , 62 ) 
oscillating plane. The mapping of the vectorial cou¬ 
pling strength A = Vc(;o|j.q measured in a 10 x 2.5/xm^ 
horizontal area in front of the magnetic bead is repro¬ 
duced in Fig. 2 f. Inherent to the dipolar structure of 
the microbead magnetic stray held, it strongly varies in 
magnitude and in orientation ex which permits a fine 
adjustment of the vectorial coupling strength with re¬ 
spect to the eigenmodes orientations by properly nano¬ 
positioning the micromagnet. Furthermore analysis of 
the fluorescence quenching in this imaging procedure also 
permits a direct identification of the locations where the 
B held is properly aligned with the NV quantization axis 
(see Fig. 2c, 2d), which is a key requirement to ensure 
efficient optical spin state readout. Finally, the triple 
requirement of avoiding fluorescence quenching, polariz¬ 
ing the Nitrogen nuclear spin and operating with a large 
parametric coupling strength compete in determining the 
best location in space where to operate the experiment. 

Having thus fully characterized the static properties of 
the system, we now investigate the qubit dynamics in 
presence of coherent mechanical motion, generated by a 
modulated piezoelectric driving force. We first restrict 
our analysis to one single mechanical mode (m = 2) by 
positioning the gradient source at a location allowing a 
large parametric coupling strength along the 02 orien¬ 
tation and tuning the external drive frequency Od/27r 
to the resonance of the second eigenmode (6.29MHz). 
The qubit is initialized in its ground state with laser 
illumination while the microwave power is adjusted so 
that Qr ~ r^d- Tflo subsequent Rabi evolution of the 
population (Jz{t) is shown in Fig. 3a, in absence (top) 
and in presence (bottom) of the coherent mechanical 
drive. The corresponding Fourier transforms are shown 
in Fig. 3b. While first a single decaying oscillation is ob¬ 
served, a longer lasting beating signal can be observed 
when the qubit is coherently oscillating in space, pre¬ 
senting a characteristic triplet spectrum, whose central 
component coincides with f^d^ revealing the synchroniza¬ 
tion of the spin on the oscillator dynamics. Increasing the 
oscillation amplitude ^r[f2d] from 0 to 9nm results in a 
linear increase of the triplet separation, see Fig. 3c, up to 
dujo « 27r X 5 MHz, corresponding to a 0.5 MHz/nm para¬ 
metric coupling strength (20 000 T/m equivalent mag¬ 
netic gradient), in good agreement with the measured 
coupling strength A -02 at the target position, see Fig. 2f. 
Scanning the MW power permits to illustrate the depen¬ 
dence of the triplet structure in the detuning between VLr 
and f^d: and underlines the synchronization of the qubit 
precession on the driven nanomotion naES]. 
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FIG. 3. A phononic Mollow triplet, a. Evolution of the 
spin population when the spin qubit is at rest (top) and oscil¬ 
lating in space (bottom) with a piezo driven 5 nm amplitude 
along 02. b. Magnitude of the corresponding Fourier trans¬ 
forms. A characteristic triplet structure is observed when 
the spin is oscillating at Qd ~ c. Triplet separation as a 
function of the oscillation amplitude. The fit corresponds to a 
parametric coupling strength of 0.5 MHz/nm. d. Dependence 
of the triplet structure on the Rabi frequency QR/27r detun¬ 
ing. e. The phonon field is dressing the microwave dressed 
qubit (see text). The measured Rabi oscillations represent a 
time resolved measurement of the dipole of the dressed qubit, 
whose spectrum reflects the allowed transitions between the 
phonon dressed multiplicities. The Mollow triplet appears 
when the oscillation amplitude is large enough to create a 
parametric energy modulation Sujq that exceeds the spin de¬ 
cay rate. It can be estimated to Tspin ~ 100 kHz from the 
decay time of Rabi oscillations in the vibrating case. 


These observations can be explained by a double dress¬ 
ing of the spin qubit with microwave photon and acoustic 
phonon fields as follows. The resonant interaction of the 
microwave pump field with the qubit can be described in 
the dressed states basis |±Ar), see Fig. 3e, parameterized 
by the number of excitations N shared between the qubit 
and the MW pump field 021. Under intense coherent ex¬ 
citation, the dynamics of the spin-microwave (polariton- 
like) subsystem can be formally described by a pseudo 
qubit (|±)), see SI, quantized along the microwave po¬ 
larization axis with a characteristic energy splitting of 
ClR/27r [49]. As a consequence of this rotation of per¬ 
spective in the Bloch sphere, the respective roles of the 


(Ty^az operators are consequently exchanged. Therefore, 
the phonon field parametrically coupled to the spin qubit 
(oc {a^a^)az) is now able to resonantly drive the pseudo 
qubit, if the resonance condition CIr ^ is met. This 
second interaction can similarly be described by a sec¬ 
ond dressing of the pseudo-qubit by the phonon field. 
This gives rise to a ladder of phonon dressed states, see 
Fig. 3e, with eigenstates |±m) parameterized by the num¬ 
ber M of phononic and dressed qubit excitations. The 
energy splitting within multiplicities can be expressed 

as (see SI) AmoIIow = ((^^d - 

which simplifies to Jci;o[^d]/2 = |(5r[Ud] -^1/2 when the 
phononic dressing field frequency Ud/27r is resonant with 
the dressed qubit energy splitting VLr/2t:. As a conse¬ 
quence the spectrum of Rabi oscillations is peaked at 
frequencies corresponding to the allowed transitions for 
the (Jz operator (see Fig. 3e). 

Measuring the temporal evolution of the spin qubit 
population {(7z{t)) indeed permits to record the tem¬ 
poral evolution of the dipole of the MW-dressed qubit 
((+|cr|—)) [49], see SI. This dipole governs the dressed 
qubit emission (in analogy with the atomic case), 
whose spectrum (see Fig. 3b) reflects the cascade among 
phononic dressed states. This situation is precisely the 
one permitting the observation of a Mollow triplet in 
QED when the atomic fluorescence spectrum under in¬ 
tense illumination was measured [25] . An important dis¬ 
tinction is that here the time resolved evolution of the 
’’atomic” dipole (the dressed qubit) is accessible. 

In order to fully explore the vectorial character of 
the parametric interaction, we now sweep the drive tone 
across both mechanical eigenfrequencies. This permits 
moving the qubit in both directions in the (01,62) os¬ 
cillating plane. For each drive frequency, the microwave 
power is adjusted to reach the resonant condition VLr « 
Ud- The measured Mollow triplet spectra are acquired 
and shown in Fig. 4a. The central component of the 
triplets is locked onto the drive frequency Od? while the 
splitting of the Mollow triplet presents two maxima, cor¬ 
responding to the response of each eigenmode. In or¬ 
der to properly understand the observed signature, it is 
necessary to precisely determine the spatial trajectories 
followed by the moving spin qubit. To do so, an opti¬ 
cal measurement similar to the one shown in Fig. Ic per¬ 
mits establishing the local orientation ep and magnitude 
5F of the electrostatic driving force field [24] and deter¬ 
mining the driven trajectories (5r(t) = 
using eq. 0 . The slight spectral overlap between the 
eigenmodes leads to elliptical trajectories (see Fig. 4b) 
which explore the oscillation plane and the magnetic field 
gradient over nano metric distances. Finally, the Mollow 
triplet’s full splitting can be adjusted with: 


(5ci;o[^] 


^ Xn.[n]{6F[n]-e^){e^-X) 

m=l,2 


(3) 
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FIG. 4. Bimodal Mollow triplet, a. Mollow triplet spec¬ 
tra obtained while sweeping the frequency f^d/27r of the elec¬ 
trostatic driving force across the two eigenfrequencies, reveal¬ 
ing the relative contributions of each eigenmodes (ei,e 2 ) to 
the Mollow triplet structure. The solid lines are fco [f^d]/2, f^d 
and Qd i (5ct;o[f^d]/2 derived from expression §, using the 
determined eigenmode orientations ei, 2 , force orientation ep 
and the vectorial parametric coupling strength Vc^o measured 
at the position marked by O in Fig. 2f. These quantities are 
reported in b where the spatial trajectories followed by the 
spin qubit are represented for varying drive frequencies. Also 
shown are the qubit iso-energy lines and the characteristic 
length StmoHow (see text), c. The normal projections of the 
elliptical trajectories on the parametric coupling vector ex 
determine the Mollow triplet separation. For illustration, the 
contribution of three specihc trajectories are highlighted in 
panels b and c (z to in). 


Both the deduced magnitude (|A|/27r = 0.5 MHz/nm) 
and the orientation (ey reported in Fig. 4b) of the vec¬ 
torial coupling constant A are in good agreement with 
the static measurements described above, at the position 
marked in Fig. 2f. Geometrically, the magnitude of the 
parametric coupling strength corresponds to the length 
of the projection of the ellipses on the ey axis. Pursuing 
this geometrical approach, it is possible to introduce a 
characteristic length, StmoWow = Tspin/IV cjqI ~ 200pm 
here (reported in Fig. 4b), which represents the minimum 
oscillation amplitude along ey necessary to resolve the 
phononic Mollow triplet. It is interesting to point out 
that this quantity is comparable to the spatial spreading 
of the nanowire Brownian motion of ~ 52 pm, responsi¬ 
ble for an equivalent incoherent parametric modulation 
of = AAxth ~ 27r X 25 kHz, which could alter the 
Mollow triplet structure [50] in larger magnetic field gra¬ 
dients. Understanding the coherence properties of the 
dressed qubit and the contribution of Brownian motion 
will be the subject of future investigation. 

Conclusion- We have demonstrated the observation of 
a phononic Mollow triplet in a spin-nanomechanical hy¬ 
brid system, reproducing with phonons and a spin qubit 
one of the founding signatures of QED based on photons 
and atoms. The observed signatures also demonstrate 
the synchronization of the spin precession onto the me¬ 


chanical oscillation frequency. This opens the road to¬ 
wards new detection strategies for observing dynamical 
spin-dependent forces SF^pinit) = In partic¬ 

ular the phononic Mollow triplet structure should also be 
imprinted on the oscillator dynamics through the spec¬ 
trum of (Jz[Ct]^ as mechanical sidebands. Pushing further 
the analogy with QED [49] it should also provide a direct 
mechanism for creating single phonon sources through 
cascaded phonon emission within the dressed state lad¬ 
der [30] . 
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